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Abstract 

Inclusive production of a quarkonium in hadron collisions at low transverse momentum can 
be used to extract various Transverse-Momentum-Dependent(TMD) gluon distributions of hadrons, 
provided the TMD factorization for the process holds. The factorization involving unpolarized TMD 
gluon distributions of unpolarized hadrons has been examined with on-shell gluons at one-loop level. In 
this work we study the factorization at one-loop level with diagram approach in the most general case, 
where all TMD gluon distributions at leading twist are involved. We find that the factorization holds 
and the perturbative effects are represented by one perturbative coefficient. Since the initial gluons 
from hadrons are off-shell in general, there exists the so-called super-leading region found recently. 
We find that the contributions from this region can come from individual diagrams at one-loop level, 
but they are cancelled in the sum. Our factorized result for the differential cross-section is explicitly 
gauge-invariant. 


1. Introduction 

Theoretical predictions of the inclusive quarkonium production at large transverse momentum in 
hadron collisions can be made by using the standard collinear factorization of QCD. Using the predictions 
one can extract from experimental results the gluon distribution functions of initial hadrons. These 
distribution functions are important for making predictions of other processes and for providing the 
information about inner structure of initial hadrons. However the information is limited because the 
extracted gluon distributions are one-dimensional. It is possible to extract three-dimensional gluon 
distributions, called as Transverse-Momentum-Dependent(TMD) gluon distributions, by using processes 
involving small transverse momenta. For this purpose, one needs to establish TMD factorizations. 

In TMD factorization nonperturbative- and perturbative effects in a process are consistently sepa¬ 
rated. TMD parton distributions are dehned with QCD operators and represent the separated nonpertur¬ 
bative effects. TMD quark distributions can be extracted from processes of Drell-Yan and Semi-Inclusive 
DIS. For these two processes, TMD factorizations have been established in mm and in 13 ll- It is 
noted that TMD factorization has been first established for inclusive e"’“e“-annihilations into two nearly 
back-to-back hadrons in the seminal work in [5], where only TMD parton fragmentation functions are 
involved. For extracting TMD gluon distributions it is suggested to use processes in hadron collisions 
like Higgs-production[6l [7], quarkonium production [8], two-photon production [9], the production of a 
quarkonium combined with a photon [TO] and double-quarkonium production[Tl]. In [121113] the TMD 
distribution of linearly polarized gluons inside a nucleus and its phenomenology have been studied. 
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TMD factorization for the suggested processes has been derived only at tree-level, except the process 
of Higgs- and quarkonium production. At one-loop level, TMD factorization has been examined for 
Higgs-production in and for the production of in M- In these studies one takes the incoming 
gluons from the initial hadrons as on-shell and spin-averaged. The corresponding perturbative coefficient 
is determined at one-loop. Since only on-shell gluons with spin averaged are considered, one examines in 
fact the TMD factorization for the contribution from the scattering of unpolarized gluons coming from 
unpolarized hadrons. By taking the transverse momenta of incoming gluons into account, the gluon as 
a parton from an unpolarized hadron can be linearly polarized according to m- Certainly one can use 
on-shell-gluon scattering as in the studies of PE] to study the factorization of the contribution from 
linearly polarized gluons at one-loop. But this is difficult because one needs to study the on-shell gluon 
scattering at three-loop, as discussed in [14]. The factorization for the processes with polarized hadrons 
has not been examined beyond tree-level. 

We will use the subtractive approach to examine TMD factorization for the inclusive production of a 
^S'o-quarkonium rjQ, which is r]c or rji,. The approach is based on diagram expansion at hadron level and 
explained in [HI |T3 E] . In this approach the nonperturbative effects can be systematically subtracted in 
terms of diagrams. Using the approach TMD factorization for Drell-Yan processes has been examined 1 19]. 
where TMD quark distributions are involved. In the case with TMD gluon distributions the situation 
is rather subtle, e.g., there can be super-leading-power contributions related to gluons as found in m, 
and TMD factorization can be violated in certain processes as shown in [snumiii]. We notice here that 
gluons from initial hadrons participating in a hard scattering are in general off-shell, i.e., the gluons have 
momenta slightly off-shell and unphysical polarizations. This brings up complications in examining the 
factorization and for obtaining gauge invariant results. The complications do not appear in [6[[T4]. In 
the subtractive approach the complications can be correctly addressed. 

In our work, we take the initial hadrons as arbitrarily polarized. We can show that at one-loop level 
the factorization holds and is explicitly gauge invariant. The contributions from the super-leading region 
do not appear at tree-level in our case, but they do appear at one-loop level. The nonzero contributions 
come from two different sets of diagrams. They are canceled in the sum. Our result shows that there is 
only one perturbative coefficient for all various contributions involving different TMD gluon distributions. 
The coefficient is determined at one-loop. 

The above discussions are mainly relevant to the initial hadrons. For the nonperturbative effects 
related to the quarkonium in the final state, one can employ the factorization with nonrelativistic 
QCD(NRQCD) proposed in [23] . A quarkonium can be taken as a bound state of a heavy quark Q 
and a heavy anti-quark Q. The heavy quark or heavy anti-quark moves with a small velocity v in the 
rest frame of the quarkonium. For the production of r]Q, the production rate at the leading order of v 
can be written as a product of the production rate of a QQ pair in color- and spin singlet with a NRQCD 
matrix element. The NRQCD matrix element characterizes the transmission of the produced QQ into 
r]Q. In this work we will also add the correction at the next-to-leading order of v in NRQCD factoriza¬ 
tion at tree-level, since the correction and the one-loop correction can be of the same importance. We 
notice here that in NRQCD factorization for a P-wave quarkonium one needs to consider not only the 
color-singlet QQ pair, but also the color-octet QQ pair [23]. Taking the color-octet QQ pair into account, 
TMD factorization for the production of a P-wave quarkonium is violated|24j. It is noted that NRQCD 
factorization with the color-octet QQ pair is also violated at two-loop level. But this violation can be 
avoided by adding gauge links in NRQCD color-octet matrix elements [25 j. 

Our work is organized as in the following: In Sect.2. we study TMD factorization at tree-level. 
Notations are introduced. In Sect.3. we study TMD factorization at one-loop level. We will show that 
one needs to introduce a soft factor to complete TMD factorization. The mentioned complications from 
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unphysical incoming gluons will be explained in detail. In Sect.4. we will give our main result of TMD 
factorization for the differential cross-section with arbitrary hadrons in the initial state. The differential 
cross-section is given in detail in the case that one of the initial hadrons is unpolarized and another one 
is of spin-1/2. Sect.5. is our summary. 

2. TMD Factorization at Tree-Level 

In this section we first introduce notations and TMD gluon density matrix or distributions in Sub¬ 
section 2.1. In Subsection 2.2. we consider the tree-level contribution with one-gluon exchange. The 
contribution from two-gluon exchange is studied in Subsection 2.3. 

2.1. Notations and TMD Glnon Distribntions 

For our purpose it is convenient to use the light-cone coordinate system, in which a vector is 
expressed as = (a"'', a“, a_L) = ((a*^ -|- a^)/\/2, (a° — a?)/\/2,a^,a^) and = (a^)^ -|- (a^)^. We 
introduce two light-cone vectors: n and 1. With these light-cone vectors one can build two tensors in the 
transverse space. The vectors and tensors are: 

^ 1 ^^ = (0,1,0,0), Z^ = (1,0,0,0), ef = (1) 

with = —e^l = 1. 

We consider the process 

hA{PA) + hsiPB) VqIq) T 

where rjQ stands for or ry^. It is a ^S'o-quarkonium consisting of a heavy QQ pair with Q 
The momenta in Eq.Q are given by: 

■P4 ~ ~ (0,T//,0,0), = {q^,q~,q±) = {xPX,yPE,Q±), (3) 

where we have neglected the masses of hadrons, i.e., P^ ~ O and ~ 0. The mass of the quarkonium 
is and the invariant mass of the initial hadrons is s ~ 2P^Pg. At the leading order of the velocity 
expansion one has ~ 2Mq. We will also use the notation q^ = = M^. We are interested in 

the kinematical region of cy^/Q <C 1. In this region, one can establish TMD factorization to express the 
differential cross-section in terms of TMD gluon distributions. 

TMD gluon distributions are defined with QCD operators. To define them for Ha we introduce the 
gauge link along the direction ,0, 0) with «+ <C 

Puiz) = P^—igs J dXu ■ A{Xu + z)^ , (4) 

where the gluon field is in the adjoint representation. The gluonic density matrix of Ha is defined as: 

with = (0 ,^“,£l) and = {xP^,0,k±). is the field strength tensor. The definition is gauge- 
invariant. The gauge links in Eq.([5]) are taken off light-cone. This is to avoid light-cone singularities if 
we take u'^ = 0. In this work we will use Eeynman gauge. In this non-singular gauge fields at infinite 
space-time are zero. If one works with a singular gauge, one needs to implement gauge links along the 
transverse direction at = —00 to make the definition gauge-invariant [26l [27]. 


(2) 

= c or Q = b. 
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The density matrix is defined for hadrons with an arbitrary spin. In general one can decompose the 
density matrix with scalar functions. These functions are TMD gluon distributions. Here we briefly 
discuss the case for Ha with spin-1/2. We assume that the spin of Ha is described by the helicity Sl and 
the transverse spin s^. At leading power or leading twist the indices /i and v of the density matrix are 
transverse. The decomposition reads jT5]: 


r7(^) = + 


2M\ 

1 


kx) + ]^Sl 


- ie^ AGL(x,k±) 


+ 


1 


2Ma . 


— k± ■ sxGTix, k±) + ie^kx ■ sxAGt{x, kx) 


i (k^^s'l + s^^k'2) i^^AHT{x,k±) - ^^AH;^{x,k±) 


• s±AHTix,k±) 


( 6 ) 


with = e^k±p and = e^sx^. The TMD gluon distribution fg corresponds to the standard gluon 
distribution function of an unpolarized hadron if we integrate out the transverse momentum formally. 
By taking the transverse momentum into account, the gluon as a parton in an unpolarized hadron can 
be linearly polarized indicated by H-^ in Eq.(|S|). 

Similarly, one can define the gluon density matrix T^^ (k) or TMD gluon distributions of where 
one uses instead of the gauge link along the direction v with v~ v~^. In the following subsections 
and Sect.3 we will set u = n and v = I for convenience to study the factorization. It is well-known that 
there will be light-cone singularities in TMD gluon density matrices with the setting. But this will not 
affect our analysis, one can always make in each step the substitution n —>■ u and I —>■ v for subtractions. 
We will make the substitution in our final result. 




Figure 1: (a): Tree-level diagram of one-gluon exchange, (b): The bubbles in the middle of Fig.la stand 
for the amplitude of g*g* — gq. 


2.2. TMD Factorization of One-Gluon-Exchange Contributions 

At tree-level the contribution to the differential cross-section of Fq. ([2]) are from diagrams represented 
by Fig.la. It can be written as: 

^ J d'^kAd'^ksSiq'^-4m^Q){2TT)‘^6'^{kA + kB-q)M‘^^f,{kA,kB)(^Mf^{kA,kB))'^ 


4 





























dH 

{2^Y 




dYrj 

w? 


e*’?-^«(/iB|^'^’^(0)A^’^(r/)|/iB) 


(7) 


In Eq.([7]), the correlation function of gluon fields in the first [• • •] of the second line is represented by the 
lower bubble in Fig.la, the correlation function in the second [• • •] is represented by the upper bubble. 
The bubble in the left-middle part of Fig.la is the amplitude given by Fig.lb, i.e., M^^{kA,kB) is 
the amplitude of g*{kA)g*{ks) Vq- At amplitude level, there is only one gluon exchanged between 
bubbles. There can be more exchanged gluons. The case of two-gluon exchange will be studied in the 
next subsection. 

We will take the leading order of the small velocity expansion in NR.QCD[23j. At the leading order, 
the quark Q and Q carry the same momentum q/2 in Fig.la. The broken line there is the cut. The 
cut cutting the lines of the heavy quark pair means also that we take the projection for the pair into 
the color-singlet ^5o-state. The projection is standard and can be found, e.g., in [29]. Under the above 
approximation, one easily finds the amplitude given in Fig.lb: 




/l(0,0) 


r{Q) 

mQ^2NcmQ' 


( 8 ) 


where we introduced the form factor fL for two-gluon fusion into r]Q. The two gluons are in general 
off-shell. The explicit result of /^(0,0) for on-shell gluons is given. Here we use the wave function at the 
origin for the projection. In the final result it will be replaced with NRQCD matrix element. 

In the kinematical region under our consideration the produced r/g carries a small transverse momen¬ 
tum q± ~ AQ with A <C 1. The leading power contributions arise when upper- and lower bubbles are 
jet-like functions. The power counting for the momenta carried by the gluons leaving the bubbles are: 


A:^~(1,A2,A,A), A:^~(A2,1,A,A). 


(9) 


The power counting for the gauge vector field in the correlation functions represented by the lower- 
or upper bubble is the same as for the momenta given in the above, respectively. Taking the tree-level 
result in Fq.(|8|) we have the leading order contribution in A for the combination appearing in Fq.(|7I): 

Ml^^ikA, kB)A^’^m^’^iv) = <?'/ l ( 0 , {k+Al^iO - kA±aA^’^iO) 

• {kBA^iv) - kB±pA^’-{v)) + 0{x% ( 10 ) 

the leading order is at A^. In the expansion one also needs to expand fL{k\,k‘^) in A. This gives the 
factor /l( 0, 0) in the above. From the power counting, the so-called super-leading region can appear here 
in the amplitude if there is a contribution proportional to A°'’~^A^'~ (rj). Such a contribution can be at 
order of A^ or A^. Usually, if we take the gluons from as on-shell, as the explicit calculation in |14j . 
one will not meet the super-leading regions, because the on-shell gluons are always transversely polarized. 
In Fq. (|inh the contribution proportional to A°''~^{^)A^'~(rj) is combined with two transverse momenta. 
Therefore, the super-leading region gives no contribution here. We also note that the introduced form 
factor in general is not gauge-invariant, but /l(0, 0) is gauge-invariant obviously. 

Since we are only interested in the leading order of A, the —components of gauge fields in the 
correlation function of Ha can always be neglected. The —components of momenta carried by these 
fields can always be neglected except in the correlation function. The latter results in that one can 
perform the integration over those —components of momenta and the corresponding -|—components 
of the space-time coordinate vectors in the correlation function trivially. Therefore, we introduce the 
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notations which will be used through this work: For any gauge field in the correlation function 

of Ha the space-time coordinate vector is = (0,The momentum k carried by this field is 
= {k~^,0,k±). Similarly, we also use these notations for the correlation function of hs in which the 
role played by )-components is exchanged by that of — (+)-components, respectively. 

With the above result the one-gluon contribution from Fig.la at the leading power of A is: 


da 


7r6{xys-Q‘^) 4 2 


s 

d^i 

(27r)^ 


9.VI(0,0) J d^kA±d^kB±6^{kA± + kB± - 


{hA I (0) {i)\hA) 


d^ff 

(27r)^ 




,( 11 ) 


with k^ = 0 and k^ = 0. Approximately one may write the two correlation functions in the second line 
as the gluon density matrix F^ and F^, respectively. However, in these correlation functions, given 
by = d^G'^ — d'^G^ is not exactly the field strength tensor operator. To obtain it, one needs to 
consider the contributions from the exchange of two- or more gluons 



Figure 2: (a) The diagram with two-gluon change. The middle bubble in the left part is the sum of 
diagrams given in (b) and (c). (b) The tree amplitude for g*g*g* gq, in which all gluons are attached 
to the quark line, (c) The tree amplitude for g*g*g* gq, in which a three-gluon vertex is involved. 


2.3. TMD Factorization of Two-Gluon-Exchange Contributions 

Now we consider the contribution in which two gluons come from Ha as in those diagrams given in 
Fig.2a. The contribution after neglecting the unimportant components of momenta and with the notation 
introduced before Ea. IfTT]) is: 


da 

d^q 


^ J d^kid?k2d^kBd{q^ - AmQ)5‘^{ki + k2 + kB- k2, kB) 


{Mf,{kl + k2,kB))^ 


d?g 

(27r)^ 


d^-^^{hB\A‘^'''{<d)A^’>^{g)\hB) 


( 12 ) 


The amplitude ^ 2 , ^s) ol contributions from Fig.2b and Fig.2c. It is the amplitude 

for three-gluon fusion into gq with the exclusion of those diagrams in which the gluon 1 and gluon 2 
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combine into one gluon with the three-gluon vertex. Those diagrams are in fact included in the lower 
bubble. A factor 1/2 has to be implemented to avoid double-counting. 

The three-gluon amplitude given by Fig.2b and Fig.2c will be used extensively in our work. Hence, 
we give the detailed result for the expansion with the power counting in A. In the expansion one should 
keep the ie factors in propagators, because these factors have a physical meaning as we will see. The 
calculation of diagrams are straightforward, although it is tedious. But after the expansion the result 
takes a rather simple form for combinations with gauge fields appearing in Eq. (|12p . We have from Fig.2b 
and Fig.2c the following: 




2b ■ 


d iid ^2d 

(27r)9 






in) 


— {B~ — term). 


— l£ 


^2 — is. 


I 

MZatih,k2,kB) J 


2c. 




G“i’+“(6)A“2’+(6)- 


1 1 

-f 


kt — is k't -\- i 


IS 


‘ ^ -i 


- + 7^— G“”+“(6)A«'+({,) 


ki -\- is k 2 — is J 


+ (B — term). 


(13) 


where the {B~ — term) is proportional to A^~{ri). It is exactly cancelled in the sum of Fig.2b and 
Fig.2c. The above results are symmetric in exchange of the gluon 1 and 2. Taking the sum we obtain 
the two-gluon-exchange contribution from Fig.2a at the leading order of A: 


da 


—^5f/L(0,0) J d^kid^kBS^{ki + kB - q)e±ati 


I Qd,-u I 


( 27 r 




( 27 r 




- i I d3A:2^^^e*^2-(6-6)_l_Al«2,+ (|2)G“i’+"(ei) 


'(2vr) 


kt + 


IS 


\hA)\. (14) 


We note here that the Tie’s in the denominators come from different places. The factor +is comes from 
the gluon propagators in Fig. 2c and indicates that the interactions through the exchange of the collinear 
gluon are of the initial-state, while the factor —is comes from the quark propagators in Fig.2b and Fig.2c 
and indicates final-state interactions. In the sum the effects of final-state interactions are completely 
canceled and there is no term with the eikonal propagator l/(A:j ^2 ~ i^)- The sum only contains those 
eikonal propagators l/(A :/'2 -fie) representing initial-state interactions. This result is important here for 
the factorization. If the effects of initial-state- and final-state interactions exist simultaneously, they can 
be potential sources for breaking the factorization, as discussed in [20l |2T], or one needs to introduce 
non-universal TMD gluon distributions for the factorization [28j . Then, the prediction power is lost. 
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Now we consider the product in the gluon density matrix defined in Eq.([^: 

I + 0(5?). (15) 

Comparing this expression with Eq. dldl) . one can realize that the first term in Eq. (jl4ji is the missing part 
for the field strength tensor in the one-gluon-exchange contribution, and the second term forms 
a part of the gauge link. One may consider the exchange of more gluons to obtain full gauge links. 
Therefore, we can write the result of TMD factorization at tree-level as: 

^ ^ gs/L(0>0) J d^kAj_(fkB±6‘^ikA± + kB±-q±)€±af,€±p^ 

r^/{kA)r'^B^{kB) (^1 + 0{as) + 0(A)) (16) 

with = {xP^,0,kA±) and k^ = {0,yPB In the above, the gluon density matrices are defined 

gauge-invariantly. /l( 0, 0) is gauge invariant as discussed before. Hence, the result in Eq. lfTH]) is gauge 
invariant. This result can be still represented by Fig.la, where the lower bubble represents r^“, the 
upper bubble represents T^^. The obtained result in Eq. ifTHl) will be corrected beyond tree-level. 


3. TMD Factorization at One-Loop-Level 

One-loop correction comes from diagrams in which one has one-gluon-exchange in the middle part 
of Fig.l. The one-loop correction can be divided into two parts. One part is the real part, in which an 
additional gluon is exchanged in the middle part of Fig.la crossing the cut. Another part is the virtual 
part, in which the exchanged gluon does not cross the cut. We will consider the two parts in the following 
subsections separately. 





Figure 3: Diagrams of the real part of one-loop correction. 


3.1. The Real Part 

The real part is given by the diagrams given in Fig 131 where the bubbles in the middle represent the 
three-gluon- and two-gluon amplitudes explained in the last section. For the bubbles attached with three 



































gluons, we always identify the two gluons coming from the below to the bubbles in the middle part are 
the gluon 1 and the gluon 2, as specified in Fig.2c. With this identification, The total contribution of the 
real part is: 



(17) 


In the kinematical region of q± ~ Q\ with A <C 1 the contributions at the leading order of A come from 
the region of the gluon momentum which is collinear to or to Pb, and soft, with the standard power 
counting. We first consider the case that the exchanged gluon is collinear to Pa in each diagram. 

The leading contribution from FigjSti can be easily derived by using the leading result of the two- 
gluon amplitude in Ea. dlOp and by taking the pattern of the momentum k of the exchanged gluon as 
~ (1, A^, A, A). The leading contribution from FiglSti is: 


^ ^ gf/L(0>0) J d^kid^kB^^6^{ki + kB-q-k){-2T:5{k^))T)^{kB) 

kBS^'^{{ki - - nai(A;i - ky/^ - kyg""- n^yki - ky^ 

+ kif^gpP + {-2ki + kYg^P^ + {ki - 2kf gP^^'^ 

-0^p^-’‘yhA\A'pyo)Aiyi)\hA) + ---. (is) 


In the above we have in fact included those contributions from exchange of all possible gluons between 
the upper bubble and the bubbles in the middle of Figj3^. Therefore we have now in the above Vy 
representing the upper bubble. This will be implicitly implied in the whole analysis. The • • • stand for 
power-suppressed contributions which are neglected. 

The power of g^ indicates that the contribution in Eg. dlSp can be a part of the 0(as)-correction 
in the tree-level factorization in Eo. dlhll . However, this needs to be examined. We note here that the 
Vy represented by the lower bubble in Eig.la is a jet-like correlation function, which is the sum of all 
diagrams in which parton lines carry the momenta collinear to Pa- Because the exchanged gluon in EiglS^ 
is collinear to Pa , the contribution from the exchanged gluon can be already included in the lower bubble 
Vy entirely or partly. Therefore, the correct contribution to the 0(as)-correction is only obtained after 
subtracting the corresponding contribution from . If one simply takes the contribution from Eig.3 as 
the 0(as)-correction in Eq. ffTHp without the subtraction, a double-counting happens. 



Figure 4: Diagrams of the gluon TMD distribution for the subtraction 
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Beyond tree-level, the TMD gluon density matrix receives the contributions from diagrams in FigjU 
In these diagrams, the double lines represent the gauge links. The contribution from FigH^ is; 

rr{t4)[^ = ^/ -h+ m}{kA ■ ns- - {h - k}V‘) 

+ (2ti - + (2k - 

{-i){kA ■ ng^- {ki - kYnp^) 

+ (-2fci + *)'’/'’> + (ii - 2kf!fY 

I (O)A^ (0\hA). (19) 


Comparing Ea. (|18l) with Ea. (ll9l) . we find that the contribution from Figj3^ takes a factorized form: 


^ 1 ^ = 0 ) J d^kA±d^kB±S\kA± + kB± - q±) 

e±.f2e±0.rY{kB)(r/^A"{kA)\^. ( 20 ) 

This result indicates that the contribution from Figl3^ is already included in the gluon density matrix 
F^“ in Eq. ffTGl) . Therefore, it should be subtracted from the C>(Q; 5 )-correction. This leads to that the 
contribution from Figl3^ will not contribute to the 0(Q;s)-correction. 

Now we consider the contribution from EiglSb. Taking the exchanged gluon as collinear to Pa and the 
leading result for the two-gluon amplitude represented by the right bubble in the middle, the contribution 
reads: 


da 


^ f d^kid^kBj^^6^{ki + kB - q - k){-2Tr6{k‘^)) 
s J ( 2 vrj^ 

d^( 




e±i3iuiki - kY - ni3^€±au{ki - kYj_ 


jfi-kB/u\nd-u 


(27r) 


{hB\G^^-Y0)A'^'Yv)\hB) 


{-ki - + ^2k - kifg'^^^^ + {2ki - kY^g^^^ 


{ki — kY — i£ 


( 21 ) 


In the above there are still some contributions at higher order of A. These contributions need to be 
separated and neglected as those represented by • • •. We note here that among diagrams in FiglSl the 
contribution from FiglSj) is the most difficult to analyze. 

To find the leading contribution we can use the leading order result of the three-gluon amplitude in 
Ea. ()13ll . We expand the relevant combination in A: 


MZatikl,-k,kB) 




Q 


(2vr) 

d^^d^fj 


m 


(2vr) 


^i^.kk+in-kB^YGk- (fj) 


(A:+ - k+)AY{i) - ki±^A^^’+{i) 


+ 


^^^Al“i’+(|) (^gBa,Ykt - k^) + + ^a.(0(A)) + 0(A2)|, 


( 22 ) 
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where the gluon field is from the correlation function of /is, and is from that of h^- 

The combination is a vector with the index 02 which will be contracted with the remaining terms in 
Ea. ()21l) . This index is carried by the exchanged gluon. It is noted that the power for different 02 is 
different. The terms proportional to in {• • •} are at order of A, while the term with 02 =T is at 
order of A^. At first look one may neglect the terms proportional to But one can not neglect them, 
because all terms contracted with the remaining terms in Eq. (|21l) will give the leading order contribution 
to the differential cross-section, except the terms proportional to Zqj. 

It will be lengthy to give the full result for the contribution from Figl3)3. However, one can rather 
easily find the factorized form of the contribution. We consider the contributions from FigHb and FigHt 
to the TMD gluon density matrix: 


r^AikA) 


H) 




J -ki + k){i){kA ■ - k'in^) 

- r—{9snpr''^^)i-i){kA ■ ng^p - {ki - kYnp^) 

+ (- 2*1 + krs^l>' + {h - 2kfgi>Y 

1 r d^kd'^kA 

xP+ J (27r)4 + k){i){kA -ng^ p + k'^Up) 

- . i9sn‘^f^'"Yi-i)ikA ■ ng^^ - {ki - kYnpY 

‘ + ti)*/-’ + (-2*1 + J:)'/* + (ki - 2k'f!fY 

I ^e‘k'‘'{hAlAfiO)A'S(()lhA)- (23) 


With these expressions and result in Eq. (j22|) . we find that the contribution from FiglSb takes the factorized 
form: 


da 

d'^q 


L 


xyndi^xys^ Q ) g4y2^Q^Q^ J (YkA±d‘^kB±d‘^{kA± + ksx - q±) 


e±apexyurYikB)(r^A"ikA)\^ 



(24) 


In the above the contribution factorized with FigHb is from these terms proportional to l/{k~^ — is) in 
the {• • •} in Eq. (l22l) . while the contribution factorized with FigHJ: is from these terms proportional to 
l/{kf +ie). Again, the contribution from Figl3)3 is totally subtracted. It does not contribute to the 
C)(Q:s)-correction. 

The remaining diagram which needs to be studied is FigJSh. The contribution involves the three- 
gluon amplitudes only. The leading order result can be obtained in a straightforward way by using the 
result in Ea. (|22n . The leading contribution comes from the term with the index 0:2 =T. We obtain the 
contribution at the leading order of A: 


da 

d'^q 


L 


^ ^ at fliO, 0 ) J dr'kid^kB-^^d'^iki + kB-q- k){-2-K5{kYYxafiexf}ukj^ 
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OLp pa\a2b pc\a2d / {ki - fc)+ 

9^} ! 


2 r 


d?i 




Tj{kB) + ---{2h) 


This contribution takes a factorized form with the contribution from FigHJl to the TMD gluon density 
matrix; 

rT(fcA)|^ = ^^0^{-27rd{e))5HkA-k, + k){i){kA-n9\ + k^np)g^J^^^^r^^ 

{-i){kA • + kPnP) (;^)'/ (26) 

We have: 


da 

A 


L 


xyTrS{xys — 

2(iV| - 1) 

e±Q/.e±/3^r^''(/i;s) 


ffs/L(0>0) J d?kA±d^kBX S^{kA± + kBX - qx) 

it 


(27) 


This indicates that the contribution from FiglSt will not contribute to the 0(as)-correction. 

From the above analysis, the leading contributions from the exchange of the gluon collinear to Pa are 
already included in the gluon density matrix in the factorized form in Eg. (11611 at tree-level. Therefore, 
these leading contributions do not contribute to the 0(Q!s)-correction in Eq. (ll6ll . Performing a similar 
analysis one will also find that the leading contributions from the exchange of the gluon collinear to Pb 
are already included in the gluon density matrix . We conclude here that the leading contributions of 
the real part from the exchange of a collinear gluon will not contribute to the 0(as)-correction. They are 
correctly factorized into TMD gluon density matrices. Therefore, the possible 0(Q;s)-correction to Eq. (jl6h 
can only come from the real part subtracted with the collinear contribution, i.e., from the difference: 


da 


da 


' xy-K 5 {xys 

d^q 

R,s 

d^q 

R 

[ 2(iV2 - 


-fi'f/L(0,0) J d‘^kA±d'^kB±5'^{kA± + kB±-q±) 


e±c 




+ h.c. 


(28) 


The contribution from the above expression at the leading order of A can only come from the momentum 
region where the exchanged gluon is soft. We will study the soft-gluon contribution in the next subsection. 

Before turning to the soft-gluon contribution, it is interesting to compare the results here with the 
calculation done in [H] by replacing the initial hadrons with on-shell gluons. In |14] one can only examine 
the factorization of a part of the differential cross-section, in which only fg in Eq.([6]) is involved, i.e., 
the contribution from unpolarized gluons in unpolarized hadrons. With the employed approach here, we 
are able to examine the factorization with the entire gluon density matrices. In this approach the initial 
gluons are in general off-shell. This brings up some additional complications in comparison with the study 
in m- The complications are indicated by the fact that we need to consider additional contributions 
represented by those contributions to gluon density matrices given by FigUt and FigHji. 


3.2 Soft-Gluon Contributions and the Soft Factor 

In the kinematical region of ~ A <C 1, dominant contributions can come from the momentum 

region in which the exchanged gluon in Figl3]is soft. The soft gluon carries the momentum k at the order: 

~ (A,A,A,A). (29) 
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In analyzing the collinear contributions, e.g., the three-gluon amplitude contracted with gauge fields 
in Eq. ljl.'lj) . there are no contributions from the super-leading region in the set of diagrams given by 
Fig.2b or Fig.2c. However, we will have the contributions from the super-leading region in analyzing the 
contributions of soft gluons. We hrst consider the three-gluon contribution contracted with gauge fields 
as that in Fq. (|22]) . 

For the case of the soft gluon we have from Fig.2b and Fig.2c the contributions at the leading order 
of A: 


,3 


= -9. 




MZa%iki,-k,kB) f 

2b J 


( 27 rj- q 


(2vr) 




= 9s 




If we calculate the soft-gluon contributions to the differential cross-section by using the result from the 
set of diagrams in Fig. 2b or Fig.2c, one will find that each contribution will be at the order of A lower 
than that of the tree-level result given in Fa. (1161) . This is the contribution from the super-leading region 
as discussed in m- However, the contributions from the super-leading region are canceled in the sum of 
Fig. 2b and Fig. 2c from Fq.([30j). 

To find the contributions from the leading region, one has to expand the contribution from Fig.2b 
and Fig. 2c at the next-to-leading order of A. The contribution from each set of diagrams is very lengthy, 
but the sum takes the simple form: 


MZlt{k^,-k,kB) J 




d^id^fj 

(27r)6 ' 


,+i^-ki+i7j-kB _ 


ml{k+ -ie) 






(l + 0(A) 


(31) 


Using this result we obtain the contribution from the soft-gluon exchange of Figi3b at the leading power: 

— 9Ag^jl{{)^Q) j d^kAd^kB j -^^6'^{kA + kB-q-k±)e^„pe±^y 


rB^kB)T^:{kA) 


iV2-l 


(27 


iV 2 -l 


k+ — ie k + ie 


1 + 0(A) , (32) 


the factor in [• • •] is a part of the soft factor introduced in the below. 

We note that there are nonzero soft-gluon contributions in the subtracted collinear contributions in 
Fq. (l28p . The contribution from the exchanged gluon collinear to Pa from FiglHb given in Fq. (j24p has the 
same soft-gluon contribution as given in Fa. (1321) . After analyzing all contributions we have the difference 
in Fq. (l28]) at the leading power: 


da 


d^q 


R^s 


— ^-^9sfL{^,0) J d^kAd^kB J -^i;^S^{kA + kB - q - k_L)e±ape±p,, 
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(33) 


The difference is nonzero at the leading order of A. If we take the factorization as given in Eg. fjlOp . then 
the difference should be taken as the 0(Q;<j)-correction. However, the difference is the effect of the soft 
gluon. It should be taken as a nonperturbative effect which needs to be factorized. For this one needs to 
implement a soft factor in Ea. ()16p . 

The needed soft factor is defined as: 

with 


1 


abc rahc 


- 1 


27r5{k^)r‘>^f 


1 


1 


k~^ — ie k + ie 


+ h.c. 


S{b^) = 


N^-l 


(0|Tr \cl(b_i_, -Qo)£u{b±, -oo)Cl{6, -oo)£y(6, -oo)j |0). 


(35) 


This soft factor has 
order one has 


been introduced in the TMD factorization of Higgs production in [ 6 ]. 


S(£±) = b^(£±} + 0(a,). 


At leading 
(36) 




(h) 



Figure 5: One-loop corrections for the soft factor. The two diagrams with their conjugated diagrams give 
the one-loop correction for the soft factor with u~^ = 0 and v~ = 0 . 


For the convenience we take u~^ = 0 and = 0 here. There are corrections at one-loop. One can 
divide the corrections into a virtual- and a real part. The real part is given by Figl5^ and the virtual 
part is given by FigUja. For u"*" 7 ^ 0 and v~ 0 there are more diagrams in which one gluon is exchanged 
between gauge links along the same direction. We will come back later to the contributions from those 
diagrams. The sum of FigJSh. and its conjugated diagram reads: 


s{i±)L =-g-s 

I 5 c-l-fe.c. 


d^k 


S\i±-k^) 


2TT6{k‘^)r’^^r^^ 


iV2-l 


1 


1 


k~^ — ie k -|- ie 


+ h.c.. 


(37) 


We note that the terms in [• • •] are exactly those in [• • •] of Eq. (j33p . Now we modify the tree-level 
factorization in Ea. (|16p as: 


da 

d'^q 


xy7r6{xys^ ) g4y2^Q^Q^ J (fkA±d^kBj_d‘^ij_d'^{kA± + kB±+£±-q±) 

e±af.exfiur^/ikAmkB)Siix) (^I + 0 (a,) + 0 (A)). 


(38) 
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With this modification, one can see that the soft-gluon contribution in Ea. ()33l) is now in the soft factor. 
Combining the results in the last subsection we conclude that the real part of one-loop correction at the 
leading order of A are included in TMD gluon density matrices and the soft factor. The real part does 
not give the 0(as)-correction. Only the virtual part can give the correction. 





Figure 6: Virtual corrections. 


3.3. The Virtual Part 

The virtual part receives contributions from two sets of diagrams. One set consists of diagrams, in 
which an additional gluon is exchanged between quark lines in diagrams given in Fig.lb. The contribu¬ 
tions from this set of diagrams contain Coulomb- and I.R. divergences. These divergences are correctly 
factorized with NRQCD factorization into the NRQCD matrix element. After NRQCD factorization, the 
contributions from this set will give a part of the 0(as) -correction. Another set of diagrams consists 
of those in which an additional gluon is exchanged between a quark- and a gluon line or between gluon 
lines. This set of diagrams is given by Figl6^ and Figl6)3 with the three-gluon amplitude defined before. 
The virtual part can be written as the sum: 

d^q V d^q 


da 

G d^q 


o 


-I- /i.c.. 


(39) 


where the last term in (• • •) is the contribution from the set of diagrams where an additional gluon 
is exchanged between quark lines. In the sum of first two terms, the contribution from the gluon- 
exchange between gluon lines is double-counted. This is corrected by the third term which stands for the 
contribution from the gluon-exchange between gluon lines. As discussed, only the first three terms are 
relevant to the TMD factorization. 

The contribution from Figl6^ is given as: 


da 

A 



7r6{xys — 


■ f d^kAd^ksd'^ikA + ks - q) f ,2 ! ■ 1,2 ! ■ 

J J (27r)4 k( + iek2+ le 




(Mf,{kA,kB)f 




( 27 r 




( 40 ) 
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with ki + k 2 = kA- A factor 1/2 should be added to avoid a double counting. The contribution is 
essentially a part of one-loop correction to the introduced form factor fL{k\,k'^). It contains collinear- 
and infrared divergences. One may think that these divergences can be regularized by the off-shellnesses 
k\ and k'g. However, in order to find the leading power contribution one has to expand the form factor 
in k\ and A:^. Only the contribution with = 0 is at the leading power, because of that 

the J-function (5 ^(A:a± + ksi, — Ql) in Figl 6 ^ is at the oder of A“^. Therefore, the contribution from 
FigJSti at the leading power of A is a part of one-loop correction to the on-shell form factor /j;,(0,0). 
Similar situation also appears in TMD factorization of Drell-Yan processes and a detailed discussion 
about virtual corrections can be found in |19] . The collinear- and infrared divergences can be regularized 
with dimensional regularization. These divergences need to correctly be factorized or subtracted. 

Using the early result one can obtain the collinear contribution in which the gluon with k 2 is collinear 
to Pa- 


J (]3kAd^kBS‘^{kA +ks - q) J 


d‘^k2 —i —i 
(27r)^ kf + ie k 2 + ie 


faia2b 


le 


k\gT - 




(27r 


(^i6 ) ( 




(2vr) 


gja^a^a^ + (-A^s - kA^g^^^^ + (kA + fci<7““^) . 


The corresponding contribution to the gluon TMD distribution from Figj 6 h is: 


r^ikA) 


1 


d^k -i 


J (27r)^ k'^ + ie 


-{i){kA ■ ng"''^^ - {kA - kyrf‘A-k){kA ■ ng^^ - k^n^ 


—i 


(41) 


M2 , -i-9sf^'')^{-k - kAT^gP^ + {2kA - k)Pg^^^ + {-kA + 

{kA — ky + le 


-{dsUpf^) 


n ■ k + ie 


d^i 

(27r)^ 


e‘i-‘‘AhA\A%mi(ei\hA), 


(42) 


where k is the momentum of the gluon attached to the gauge link and is flowing into the gauge link. 
Comparing Ea. ()41l) with Eg. (1421) . one finds: 


da 

d‘^q 



xy'K6{xys — Q^) 

2(iV| - 1) 


gtf 1(0,0) 


^/{kB) (r^A^ikA) 



d^kAd?kB5'^{kA + kB - q)exap.(^xyy 


(43) 


therefore, the collinear contribution is already included in the TMD gluon density matrix. Performing 
the analysis for the case that the gluon is collinear to Pb, one obtains the similar result. We then have 
the difference 


da 


da 


'xyTrS{xys 

d^q 

V,s 

d^q 

V 

[ 2(iV2 - 


- 55 / 1 ( 0 , 0 ) J d^kAxd^kBX d^{kAX + kBX - q±) 


e±apexpJr/{kB)r'^/(kA) ^ +r^/{kA)T’y{kB) 

V Fig^ 




-|- h.c. 


(44) 
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which does not contain any collinear divergence. But, the difference contains infrared divergences from 
the soft-gluon exchange. 

For the soft-gluon exchange we consider the case that is soft, i.e., ki^ ~ (A, A, A, A). In this case, ki 
is collinear. Therefore, in analyzing the soft contribution from FiglG^, the factor 1/2 should be replaced 
with 1 to obtain the correct result. Using the result in the last subsection for the three-gluon amplitude, 
we have the soft-gluon contribution: 


^ at 0) J d^kAd^keS^ikA + ke - q)e±pf,e±fi. 


2 (iV 2 - 1 ) 

r^^{kB)rf{kA) 


9s 


- 1 


j^aia 2 CL j^aia 2 CL 


(rk2 —iii 
( 27 r)^ ^2 + ^2 — is k2 + is 


(45) 


Again, there is also a soft-gluon contribution in the gluon TMD distribution from FigjGh with the same 
factor in [• • •] in Eq. (j45ji . At the end one finds the soft-gluon contribution for the difference in Eq. (j44h 
which is similar to the case of real corrections: 

_ xyTT 6 {xys - Q‘^) 4 2, 


da 

d^q 


V,s 


2(iV2 - 1) 


9s 


- 51 / l ( 0 ) 0 ) J d^kAd^ksd'^ikA + ks - q)sxp^lSx0yV''^{k 

d^k 2 —ii i 


jaia 20 , jaia 2 CL 


iVl-l 


( 27 r)^ kl + is ^2 - is k^ + is. 


+ h.c., (46) 


This soft-gluon contribution is in fact included in the soft factor. The soft factor receives the contribution 
from FiglSh and its conjugated diagram. It is: 


S{e±,y, 



z2(o \ 9s fobc fobc 


d'^k —ii i 

{ 2 'iiY k~ — is k~^ + is 


(47) 


If we take the factorized form as in Ea. (l38p . then the soft-gluon contribution is included in the soft factor. 

Based on the results in this subsection and previous ones, we find that at the leading power of q± ~ XQ 
with A <C 1, the real part of the one-loop correction are correctly factorized into TMD gluon densities 
and the introduced soft factor. It will not contribute to the 0(as)-correction. The virtual part gives 
contributions to the correction, but the collinear- and infrared divergences are subtracted into TMD 
gluon density matrices and the soft factor, respectively. The correction is finite. 


4. The Final Result 

In the study of the previous sections, we have set u = n and v = I for convenience. This setting will 
generate light-cone singularities in the subtraction. To present the final result, we undo the setting for 
TMD gluon density matrices and the soft factor. For u~^ ^ 0 and v~ ^ 0 the TMD density matrix 
or will depend on an extra parameter (u or respectively. The soft factor contains the parameter 
p. These parameters are defined as: 


Sr 




u v'^ 


(48) 


Our final factorized result can be written as 
da 2TT6{xys — 


dxdydP'qx 


ao J d‘^kA±d‘^kB±d‘^l±6‘^{kA± + kBx + - 9±)e±Q:/^e±/3i- 


Q 2 

mCu, Cv, p)(rf^A"{kA, Cu)rB^{kB, Cv) ] S{i^,p), 


(49) 
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where the dependence on the renormalization scale ^ in each term in the last line is suppressed, 
given by 


o-Q 


{Airas)^ 

N,{N^ - l)mQ 


\m\" 


(47ras)2 

2iV2(iV2 - l)mQ 




fjo is 
(50) 


Here, we have expressed the quantity |V'(0)p with the corresponding NRQCD matrix element (©(^S'g^^)). 
The definition of the matrix element can be found in [23]. 

As mentioned in the subsection 3.2. there are more diagrams for the one-loop correction of the 
soft-factor for the case of 7 ^ 0 and v~ 7 ^ 0. In these diagrams there is one-gluon exchange between 
gauge links along the same direction. In this case the TMD gluon density matrices also receive one-loop 
contributions from one-gluon exchange between gauge links along the same direction. These contributions 
are exactly canceled in Eq. (l38h or Eq. (|49h by those from the soft factor. 

In Ea. ()49p % is the perturbative coefficient which starts at the order of Because of that all 
one-loop real corrections are subtracted into TMD gluon density matrices and the soft factor from our 
analysis in Sect. 3.2., the coefficient is determined by the virtual correction. It is determined by the 
form factor of the fusion of two on-shell gluons into rjQ after the subtraction of collinear- and infrared 
divergences with TMD gluon density matrices and the soft factor as shown in Sect. 3.3. Since there 
is only one form factor for the fusion, we have then correspondingly in Eg. (1491) only one perturbative 
coefficient. It is noted that the form factor of the fusion with two on-shell gluons and the subtraction are 
gauge-invariant, H and Ea. dl^ are hence also gauge-invariant. 

At the leading order of ^ is 1 . Beyond the leading order will depend on (u, Cv and p because 
of the subtraction. The dependences will be canceled by those of TMD gluon density matrices and the 
soft factor. TL is obtained in m- But there are several typos and errors in constant terms. We will 
give here the corrected one. For NRQCD factorization we have made the expansion in the small velocity 
V, as discussed in the section of Introduction. We have only taken the leading order v^. However, the 
correction from the next-to-leading order of v, i.e., the relativistic correction, is at the same level of the 
importance as the 0(Q!s)-correction, as discussed in [23]. Here we also include the relativistic correction. 
This correction can be extracted from the results in [29]. We have: 


4 (iP(i5W)) , a,N, 


nCu,Cv,p) = 1 - 






47r 


a2 a2 

1 2 I 1 2 1 2 


I -k 2 In 




+ 2 In 


K 




5 —-vr" 
4 


+ 0{al) + 0{v^) 


(51) 


with = 4mQ. In the second term in Eg. (15111 there is a ratio of two NRQCD matrix elements defined 
in |23j . This term is at order of The neglected correction from the expansion of the small velocity v 
is now at order of v'^. 

After giving our main result in Eg. (1491) it is worthy to discuss the problem of the so-called scheme- 
dependence in TMD factorization. The scheme-dependence arises because one can define different TMD 
parton distributions. In the case of TMD quark distributions one can define subtracted quark distri¬ 
butions to absorb the corresponding soft factor as suggested in m- One may also use the definition 
from the soft-collinear effective theory given in m- With different definitions one obtains the similar 
factorized result with different perturbative coefficients. The difference can be calculated perturbatively 
as discussed in |32ll33j . Similarly, one can also work with different definitions of TMD gluon distributions. 
In this work, we only give our result in Eo. (|49p with the unsubtracted TMD gluon distributions defined 
in Eq.([5]). Hence, the soft factor S appears explicitly. 
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The TMD gluon density matrix of Ha or hs depends on the parameter or respectively. The 
dependence is determined by Collins-Soper equation. This equation can be used to resume terms of 
the large logarithms of q±IQ in perturbative coefficient functions of collinear factorization. In this way, 
one obtains the standard Colhns-Soper-Sterman(CSS) resummation formalism in [Tj. However, in CSS 
formalism there can be certain freedom to re-define perturbative coefficient functions to make them 
process-independent or universal. This has been noticed in [M]. An application by using the freedom 
is given in [35] for Higgs production in hadron collisions. Following the work in [3l|, the impact of the 
scheme-dependence in TMD factorization on the correspondingly derived resummation formalism has 
been studied in [36|. The difference between different schemes can be determined perturbatively. It is 
noted that our scheme of TMD factorization here and hence the corresponding resummation formalism 
are referred as Ji-Ma-Yuan scheme according to [33] [36]. 

From our result one can derive the resummation formula of \^{q±/Q) for the production of rjf.fi. The 
resummation for the production of J/V' or T has been studied in [37]. In these resummations the quantum 
numbers of the produced heavy quark pair are fixed. The resummation for the production of a heavy 
quark pair in general case is studied in [38l |39] . 

In the factorization formula in Eq. (j49|) the physical effects from initial hadrons in the processes only 
appear in TMD gluon density matrices. They take different parametrization forms for different hadrons, 
e.g., different spins of hadrons. With the formula one can derive the angular distribution for given hadrons 
in the initial state. But, the results can be very lengthy. Here, we consider a realistic case in which hs 
is of spin-0 or unpolarized, Ha is of spin-1/2. The classification of the gluon TMD distributions of a 
spin-1/2 hadron is given in Eq.Q. Eor we have; 




(52) 


With the factorization formula we can derive the result of the differential cross-section in the considered 

C3;SG clSI 


dcrix,y,q±) 


2TT5{xys — 


-(XqH 


Mx, y, q±) + s± ■ q±B{x, y, 


(53) 


dxdyd?q± 

where s_j_ is the transverse-spin vector and s±u. The two coefficient functions are expressed with 

gluon TMD distributions as: 

A{x,y,q±) = ]: [ d^kA±d^kB±d^£±6‘^{kA± + kB±+i±-qj_)S{i±,p) 


1 


B{x,y,q±) = 


-1 


. 2 / d^kA±d^kB±d^£_LS^{kA± + kB±+£±-q±)S{£_L,p) 

2MAqi J 


fgix,kA±)fg{y,kB±) - [^ikA± ' ^Bx) -kA^k^AH {x,kA±)H iy,kB±) 


kAX ■ qxfgiv, kBl_)GT{x, kA±) + 

(^AHT{x,kAA - ^^AH:^{x,kA±) 

kA±■kB± ( 


1 


8M| 


2kA± ■ kBxqx ■ kB± + k%j_kAi. ■ q± 


+2 


Ml 


■ [q± ■ kA±kA± ■ kB± + kA±kB± ■ q± ) Ai7^(x, kA±) 


H-^{y,kB±) 


(54) 


At tree-level, i.e., with 7^ = 1 in Ea. ()53p . parts of the above expression has been derived before. The 
function A is given in [8] and the first term in B is derived in m- In |14j the factorization of the 
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contribution with fgfg is examined at one-loop level. In this work, we have examined the factorization 
of all contributions at one-loop with the general result given in Eg. (149 j) . 

The obtained differential cross-section does not depend on the helicity of the polarized hadron. But it 
depends on the transverse spin. This dependence will give an Single transverse-Spin Asymmetry(SSA). 
The same SSA has been studied with the twist-3 collinear factorization in [40] . In the kinematical region 
^QCD (?_L <5, both the factorizations apply, as shown for SSA in Drell-Yan processes in [H], where 

SSA is only generated by Sivers quark distribution in TMD factorization. In our case SSA is not only 
generated by the gluonic Sivers function Gt, but also by other two T-odd TMD gluons distributions. 
Therefore, the relation between the two factorizations for SSA needs to carefully examined. 

5. Summary 

We have studied TMD factorization for r/Q-production in hadron collisions at low transverse mo¬ 
menta. If the factorization holds, one can use the production process to extract TMD gluon distributions 
from experiments. We have explicitly shown that the factorization holds at one-loop level, in which all 
nonperturbative effects are factorized into TMD gluon density matrices and a soft factor defined as the 
vacuum expectation value of product of gauge links. There is only one perturbative coefficient standing 
for all perturbative effects. This coefficient is determined at one-loop level and implemented with the 
relativistic correction of r]Q. With the result here, all TMD gluon distribution functions at leading power 
can be safely extracted from experimental data. 

In general the initial gluons from the initial hadrons are off-shell, and they can be with nonphysical 
polarizations. This makes the study of the TMD factorization more complicated than that in the case 
of on-shell gluons. However, with the complication we can still show that the factorization holds and is 
gauge-invariant. It is interesting to note that at one-loop level there exist contributions from the super¬ 
leading region from different sets of diagrams. But, they are cancelled in the sum. With our factorized 
result at one-loop, it is possible to show the factorization beyond one-loop for the studied process. One 
can also use the approach employed here for examining TMD factorization of other processes involving 
TMD gluon distributions mentioned in the Introduction. 
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